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Models for Type la Supernovae (SNe la) are reviewed. It is shown that there are strong reasons 
| to believe that SNe la represent thermonuclear disruptions of C-O white dwarfs, when these 

, white dwarfs reach the Chandrasekhar limit and ignite carbon at their centers. 

' Progenitor scenarios are reviewed critically and the strengths and weaknesses of each scenario 

are explicitly presented. It is argued that single-degenerate models, in which the white dwarf 
» I . accretes from a subgiant or giant companion are currently favored. The relation of the different 

models to the use of SNe la for the determination of cosmological parameters is also discussed. 
Observational tests of the conclusions are suggested. 
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1. Introduction 

During the past year two groups have presented strong evidence that the expansion 



! °f the universe is accelerating rather than decelerating (Perlmutter et al. 1998, 1999; 

Schmidt et al. 1998; Riess et al. 1998). This surprising result comes from distance 
measurements to about fifty supernovae Type la in the redshift range z = to z = 1. The 
results are consistent with the cosmological constant (or vacuum energy) contributing to 
the total energy density about 70% of the critical density. 

This unexpected finding, as well as the use of supernovae Type la to measure the 
Hubble constant (e.g. Sandage et al. 1996; Saha et al. 1997), has focused the attention 
again on the frustrating fact that in spite of decades of research the exact nature of the 



progenitors of supernovae Type la remains unknown. Until this problem is solved, one 
+-> ■ cannot be fully confident that supernovae at higher redshifts are not somehow different 

from their low redshift counterparts. In the present review I therefore examine critically 
the question of the nature of the progenitors of supernovae Type la. Other recent reviews 
include Branch et al. (1995), Livio (1996a), Renzini (1996) and Iben (1997). 

1' 

2. Characteristics and the basic model 

The defining characteristics of supernovae Type la (SNe la) are both spectral: (i) the 
lack of lines of hydrogen, and (ii) the presence of a strong red Si II absorption feature 
(A6355 shifted to - 6100 A). 

Once defined as SNe la, the following are several of the important observational char- 
acteristics of the class which may help in the search for progenitors: 

(1) Homogeneity: Nearly 90% of all SNe la form a homogeneous class in terms of their 
spectra, light curves, and peak absolute magnitudes. The latter are given by 

Mb ^ M v ~ -19.30(±0.03) + 51og(i/ /60 km s" 1 Mpc" 1 ) (2.1) 

with a dispersion of a{M-Q) ~ cr(Mv) ~ 0.2-0.3 (Hamuy et al. 1996a; Tamman & Sandage 
1995; and see Branch 1998 for a review). 

(2) Inhomogeneity: Some differences in the spectra and light curves do exist (e.g. 
Hamuy et al. 1996b). In terms of explosion strength, SNe la can roughly be ordered as 
follows: SN 1991bg and SN 1992K represent the weakest events, followed by weak events 
like 1986G, followed by about 90% of all SNe la which are called "normals" (or "Branch 
normals"), to the stronger than normal events like SN 1991T. 
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(3) The luminosity function of SNe la declines very steeply on the bright side (e.g. 
Vaughan et al. 1995). Since selection effects cannot prevent the discovery of SNe which 
are brighter than the "normals," this implies that the normals are essentially the bright- 
est. 

(4) Near maximum light, the spectra are characterized by high velocity (8000-30,000 
km s _1 ) intermediate mass elements (O-Ca). In the late, nebular phase, the spectra 
are dominated by forbidden lines of iron (e.g. Kirshner et al. f 993; Wheeler ct al. 1995; 
Ruiz-Lapuente et al. 1995; Gomez et al. 1996; Filippenko 1997). 

(5) Fairly young populations appear to be most efficient at producing SNe la (e.g. they 
tend to be associated with spiral arms in spirals; Delia Vallc & Livio 1994; Bartunov, 
Tsvetkov & Filimonova 1994), but relatively old populations (t^4x 10 9 yr) can also 
produce them. In particular, SNe la do occur in ellipticals (e.g. Turatto, Cappcllaro & 
Bcnctti 1994). This immediately implies that SNe la are not caused by the core collapse 
of stars more massive than 8 M© . 

(6) There exist a number of correlations between different pairs of observables (see 
e.g. Branch 1998 for a review). Of these, the most frequently used in the context of 
determinations of cosmological parameters is the correlation between the absolute mag- 
nitude and the shape of the light curve. Basically, brighter SNe la decline more slowly. 
The parameter commonly used to quantify the light curve shape is Ami 5 (Phillips 1993), 
the decline in magnitudes in the B band during the first 15 days after maximum light. 
Hamuy et al. (1996a) find slopes dM B /dAm 15 = 0.78±0.17, dM Y /dAm 15 = 0.71 ±0.14, 
and dMi/dAmi5 = 0.58 ± 0.13. Using a stretch-factor s (Perlmutter et al. 1997), one 
can write M B = M B (s = 1) - a * (s - 1), with M B (s = 1) = -19.46 (e.g. Sandage et al. 
1996), and a = 1.74 (Perlmutter ct al. 1999). Sophisticated techniques for using the 
different correlations in distance determinations have been developed (e.g. Riess et al. 
1996, 1998). 

The above characteristics can be augmented with the following suggestive facts: 

(1) The energy per unit mass, l/2(~ 10 4 km s -1 ) 2 , is of the order of the one obtained 
from the conversion of carbon and oxygen to iron. 

(2) The fact that the event is explosive suggests that degeneracy may play a role. 

(3) The spectrum contains no hydrogen. 

(4) The explosions can occur with long delays, after the cessation of star formation. 
All the properties above have led to one agreed upon model: SNe la represent ther- 
monuclear disruptions of mass accreting white dwarfs. 

It is interesting that there exists a unanimous consensus on this model in spite of the 
fact that the essence of flame physics and the details of the transition from deflagration 
to detonation (in particular the density at which the transition occurs) , which are at the 
heart of the model, remain as major unsolved problems (e.g. Khokhlov, Oran & Wheeler 
1997; Woosley 1997; Reinecke, Hillcbrandt & Niemeyer 1998; and talks by Khokhlov, Ar- 
nett, and Hillebrandt presented at the Chicago meeting on Type la Supernovae: Theory 
and Cosmology, October 1998). In fact, given these uncertainties in the deflagration to 
detonation transition it is almost difficult to understand how the entire family of SNe la 
light curves can be fitted essentially with one parameter (e.g. Perlmutter et al. 1997), 
although it is possible that all SNe la explode at the same WD mass (see §4), and that 
the entire observed diversity stems from different 56 Ni masses. 

3. Why is identifying the progenitors important? 

The fact that we do not know yet what are the progenitor systems of some of the 
most dramatic explosions in the universe has become a major embarrassment and one of 
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the key unsolved problems in stellar evolution. There are several important reasons why 
identifying the progenitors has become more crucial than ever: 

(i) The use of SNe la as one of the main ways to determine the key cosmological 
parameters H , and the contribution to the energy density Qm, requires an under- 
standing of the evolution of the luminosity, and the SN rate with cosmic epoch. Both of 
these depend on the nature of the progenitors. 

(ii) Galaxy evolution depends on the radiative, kinetic energy, and nucleosynthctic 
output of SNe la (e.g. Kauffmann, White & Guiderdoni 1993). 

(iii) Due to the uncertainties that still exist in the explosion mechanism itself, a knowl- 
edge of the initial conditions and of the distribution of matter in the environment of the 
exploding star are essential for the understanding of the explosion. 

(iv) An unambiguous identification of the progenitors, coupled with observationally 
determined SNe la rates will help to place meaningful constraints on the theory of binary 
star evolution (e.g. Livio 1996b; Li & van den Heuvel 1997; Yungelson & Livio 1998; 
Hachisu, Kato & Nomoto 1999). In particular, a semi-empirical determination of the 
elusive common-envelope efficiency parameter, acE, may be possible (e.g. Iben & Livio 
1993). 

4. Refinements of the basic model 

The basic model for SNe la (that essentially all researchers in the field agree upon) is 
that of a thermonuclear disruption of an accreting white dwarf (WD). However, additional 
refinements to the model are possible on the basis of existing observational data and 
theoretical models. These refinements still do not involve the question of the progenitor 
systems. Rather, they address the question of the WD composition, and of its mass at 
the instant of explosion. 

4.1. The composition of the exploding WD 

In principle, the WD that accretes to the point of explosion could be composed of He, of 
C-O, or of O-Ne. Let us examine these possibilities one by one. 

(i) He WDs: Helium WDs have typical masses that are smaller than ~ 0.45 M Q (e.g. 
Iben & Tutukov 1985). While if accreting, these He WDs can explode following central 
He ignition at ~ 0.7 M Q , the composition of the ejected matter in this case will be that of 
He, 56 Ni and decay products (e.g. Nomoto & Sugimoto 1977; Woosley, Taam & Weaver 
1986). This is entirely inconsistent with observations (observational characteristic (4) in 
§2). Therefore, He WDs certainly do not produce the bulk of SNe la. 

(ii) O-Ne WDs: Oxygen-Neon WDs form in binaries from main sequence stars of 
~ 10 M Q , although the precise range which allows formation is somewhat uncertain 
(e.g. Iben & Tutukov 1985; Canal, Isern & Labay 1990; Dominguez, Tornambe & Isern 
1993). These systems are probably not numerous enough to constitute the main channel 
of SNe la (e.g. Livio & Truran 1992). It is also generally expected that O-Ne WDs 
that manage to accrete enough material to reach the Chandrasekhar limit will produce 
preferentially accretion-induced collapses (to form neutron stars) rather than SNe la (e.g. 
Nomoto & Kondo 1991; Gutierrez et al. 1996). I should note that the existing calculations 
have been performed for WDs of O-Ne-Mg composition, while recent calculations of the 
evolution of a 10 M Q star produce degenerate cores which are almost devoid of magnesium 
(Ritossa, Garcia-Berro & Iben 1996). Nevertheless, because of the above two points it is 
unlikely that O-Ne WDs produce the bulk of SNe la. 

(iii) C-0 WDs: Carbon-Oxygen WDs are formed in binaries from main sequence 
stars of up to ~ 10 M Q . They are therefore both relatively numerous, and they provide a 
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significant "phase space volume" (masses in the range 0.8-1.2 M Q ; accretion rates in the 
range 10 -8 -10 -6 M Q /yr) in which they are expected to produce SNe la (upon reaching 
the Chandrasekhar limit; e.g. Nomoto & Kondo 1991). Consequently, the accreting WDs 
that produce most of the SNe la are very probably of C-0 composition^. 

4.2. At what mass does the WD explode and where and in what fuel does the ignition 

take place? 

While there is virtually unanimous agreement about everything I said up to now, namely, 
that: SNe la are thermonuclear disruptions of accreting C-0 WDs, the next step in the 
refinement of the model is more controversial. Two major classes of models have been 
considered, and they suggest entirely different answers to the questions posed by the title 
of this subsection. In one class, the WD explodes upon reaching the Chandresekhar mass, 
as carbon ignites at its center. In the second, the WD explodes at a sub- Chandresekhar 
mass, as helium ignites off-center. I will now review briefly each of these classes and 
point out its strengths and weaknesses. 

4.2.1. Chandrasekhar mass carbon ignitors 

In this model, considered 'standard,' the WD accretes until it approaches the Chan- 
drasekhar mass. Carbon ignition occurs at or very near the center and the burning front 
propagates outwards. The main strengths of this model are (see e.g. Hoefhch & Khokhlov 
1996; Nugent et al. 1997 for detailed modeling): 

(1) Some 10 51 ergs of kinetic energy are deposited into the ejecta by nuclear energy. 

(2) 56 Ni decay powers the lightcurve. 

(3) The density and composition as a function of the ejection of velocity (X;(V e j)) are 
consistent with the observed spectra. 

(4) The fact that the explosion occurs at the Chandrasekhar mass explains the homo- 
geneity. 

(5) Spectra (e.g. of SNe 1994D, 1992A) can be fitted in great detail by theoretical 
models (e.g. Nugent et al. 1997). 

The main weaknesses of the Chandrasekhar mass models are: 

(1) It has proven more difficult than originally thought for WDs to accrete up to the 
Chandrasekhar mass in sufficient numbers to account for the SNe la rate. The difficulty 
is associated with mass loss episodes in nova explosions, in helium shell flashes and in 
massive winds or common envelope phases. I will return to some of these problems when 
I discuss specific progenitor models. 

(2) For initial WD masses larger than ~ 1.2 M Q , accretion-induced collapse is a more 
likely outcome than a SN la (e.g. Nomoto & Kondo 1991). 

(3) The late-time spectrum (~ 300 days), and in particular the Fe III feature at 
~ 4700 A does not agree well with Chandrasekhar mass models (Liu, Jeffrey & Schultz 
1998). 

My overall assessment of Chandrasekhar mass models is that the strengths significantly 
overweigh the weaknesses. The calculation of late-time, nebular spectra involves many 
uncertainties, and hence I do not regard weakness (3) above as fatal (although clearly 
more work will be required to explain it away). Both weaknesses (1) and (2) can be 
overcome if it can be demonstrated that SNe la statistics can be reproduced within the 
uncertainties that still plague the theoretical population synthesis models. As I will show 
in §5, this is indeed the case. 
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4.2.2. Sub-Chandrasekhar mass helium ignitors 

In these models a C-0 WD accumulates a helium layer of <~ 0.15 M with a total mass 
that is sub-Chandrasckhar. The helium ignites off-center (at the bottom of the layer), 
resulting in an event known as "Indirect Double Detonation" (IDD) or "Edge Lit Detona- 
tion" (ELD). Basically, one detonation propagates outward (through the helium), while 
an inward propagating pressure wave compresses the C-0 core which ignites off-center, 
followed by an outward detonation (e.g. Livne 1990; Livne & Glasner 1991; Woosley 
& Weaver 1994; Livne & Arnett 1995; Hoefhch & Khokhlov 1996; and Ruiz-Lapuente, 
talk presented at the Chicago meeting on Type la Supernovae: Theory and Cosmology, 
October 1998). 

The main strengths of ELD (sub-Chandrasckhar) models are: 

(1) It is easier to achieve the required statistics, since less mass needs to be accreted, 
and the WD does not need to be extremely massive (e.g. Ruiz-Lapuente, Canal & Burkert 
1997; Di Stefano et al. 1997; Yungelson & Livio 1998). 

(2) The late-time spectrum (in particular the Fe III feature at ~ 4700 A) agrees better 
with ELD models. 

(3) SNe la light curves can be reproduced adequately by ELD models (although the 
light curves rise somewhat faster than observed, due to 56 Ni heating; Hoefhch et al. 
1997). 

The main weaknesses of ELD models are: 

(1) The spectra that are produced by ELD models generally do not agree with ob- 
servations (e.g. of SN 1994D; Nugent et al. 1997). The agreement is somewhat better 
for the subluminous SNe la (e.g. SN 1991bg; Nugent et al. 1997; Ruiz-Lapuente, talk 
presented at the Chicago meeting on Supernovae, October 1998), but even there it is not 
very good. 

(2) The highest velocity ejecta have the wrong composition ( 56 Ni and He; not inter- 
mediate mass elements; also no high velocity C; e.g. Livne & Arnett 1995). This is due 
to the fact that in these models, essentially by necessity, the intermediate mass elements 
are sandwiched by Ni and He/Ni rich layers. 

(3) Since ELD models allow for a range of WD masses, and since more massive WDs 
produce brighter SNe, one might expect this model to produce a more gradual decline on 
the bright side of the luminosity function, in contradiction to the observed sharp decline 
(see §2 characteristic (3)). 

My overall assessment of the sub-Chandrasckhar mass model is that the weaknesses 
(and in particular weakness (2) which appears almost inevitable) greatly overweigh the 
strengths in terms of this being a model for the bulk of SNe la. It is still possible that 
ELDs may correctly represent some subluminous SNe la (e.g. Ruiz-Lapuente, Canal, & 
Burkert 1997). I should note that Pinto (verbal communication at the Chicago meeting on 
Supernovae) insists that his ELD models manage to overcome all of the above weaknesses 
and that they are able to produce excellent fits to both light curves and spectra. By the 
time of the writing of this review, however, I have unfortunately failed to find published 
results of these models and hence I cannot comment on them. 

4.3. The favored model 

On the basis of the above discussion the basic model can be further refined, and I tenta- 
tively conclude that: SNe la represent thermonuclear disruptions of mass accreting C-0 
white dwarfs, when these white dwarfs reach the Chandrasekhar limit and ignite carbon 
at their center s\ 
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5. The two possible scenarios 

The next step in the search for the progenitor systems of SNe la is even more contro- 
versial. Two possible scenarios have been proposed: (i) The double-degenerate scenario, 
in which two CO WDs in a binary system are brought together by the emission of grav- 
itational radiation and coalesce (Webbink 1984; Iben & Tutukov 1984). (ii) The single- 
degenerate scenario, in which a CO WD accretes hydrogen-rich or helium-rich material 
from a non-degenerate companion (Whelan & Iben 1973; Nomoto 1982). 

In the first scenario the progenitor systems are necessarily binary WD systems in which 
the total mass exceeds the Chandrasekhar mass, and which have binary periods shorter 
than about thirteen hours (to allow merger within a Hubble time). 

In the second scenario the progenitors could be systems like: (i) Recurrent novae (both 
of the type in which the WD accretes hydrogen from a giant like T CrB, RS Oph and 
of the type in which the WD accretes helium rich material from a subgiant like U Sco, 
V394 CrA, and Nova LMC 1990#2), (ii) Symbiotic Systems (in which the WD accretes 
from a low mass red giant), or (iii) persistent Supersoft X-ray Sources (in which the WD 
accretes at a high rate ^ 10~ 7 M Q /yr from a subgiant companion). 

I will now examine the strengths and weaknesses of each one of these scenarios. 

5.1. The double-degenerate scenario 

There is no question that binary white dwarf systems are an expected outcome of binary 
star evolution (e.g. Iben & Tutukov 1984; Iben & Livio 1993). Once the lighter WD 
(which has a larger radius) fills its Roche lobe, it is entirely dissipated within a few 
orbital periods, to form a massive disk around the primary (e.g. Rasio & Shapiro 1994). 
The subsequent evolution of the system depends largely on the accretion rate through 
this disk (e.g. Mochkovitch & Livio 1990; see discussion below). 
The main strengths of this scenario are the following: 

(1) The absence of hydrogen in the spectrum is naturally explained in a model which 
involves the merger of two C-0 WDs. In fact, if hydrogen is ever detected in the spectrum 
of a SN la, this would deal a fatal blow to this model. Tentative evidence for circumstellar 
Ha absorption is SN 1990M was presented by Polcaro and Viotti (1991). However, Delia 
Valle, Benetti & Panagia (1996) demonstrated convincingly that the absorption was 
caused by the parent galaxy, rather than by the SN environment. 

(2) In spite of some impressions to the contrary, many double WD systems do exist. In 
a sample of 153 field WDs and subdwarf B stars, Saffcr, Livio & Yungelson (1998) found 
18 new double-degenerate candidates. There are currently eight known systems with 
orbital periods of less than half a day. While only one of those systems (KPD 0422+5421; 
Koen, Orosz & Wade (1998)) has a total mass which within the errors could be higher 
than the Chandrasekhar mass, the sample of confirmed short-period double-degenerates 
is still smaller than the number predicted to contain a massive system. 

(3) Population synthesis calculations predict the right statistics for mergers, about 
10~ 3 yr" 1 events for populations that are ~ 10 8 yr old and 10 -4 yr -1 for populations 
that are - 10 10 yr old. 

(4) Since double WD systems were found to exist, mergers with some "interesting" 
consequences (cither a SN la or an accretion-induced collapse) appear inevitable. 

(5) The explosion or collapse is expected to occur at the Chrandrasekhar mass, which 
as I noted in §4.3, I regard as a property of the favored model. 

The main weaknesses of the double-degenerate scenario are the following: 
(1) There are strong indications that WD mergers may lead to off-center carbon igni- 
tion, accompanied by the conversion of the C-0 WD to an O-Ne-Mg composition, and 
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followed by an accretion-induced collapse rather than a SN la (e.g. Mochkovitch & Livio 
1990; Saio & Nomoto 1985, 1998; Woosley & Weaver 1986). 

(2) Galactic chemical evolution results, and in particular the behavior of the [O/Fc] 
ratio as a function of metallicity ([Fe/H]) have been claimed to be inconsistent with WD 
mergers as the mechanism for SNe la (Kobayashi et al. 1998). 

Since we are now getting to the final stages in the identification of the progenitors, 
it is important to assess critically the severity of the above weaknesses. I will therefore 
discuss now each one of them in some detail. 

5.1.1. Constraints from Galactic chemical evolution 

Supernovae Type II (SNe II) are explosions resulting from the core collapse of massive 
( 8 M©) stars. These supernovae produce relatively more oxygen and magnetism than 
iron ([O/Fe] > 0). On the other hand SNe la produce mostly iron and little oxygen. Until 
recently, the impression has been that metal poor stars ([Fe/H] < —1) have a nearly flat 
relation of [O/Fe] vs. [Fe/H], with a value of [O/Fe] - 0.45 (e.g. Nissen et al. 1994), while 
disk stars ([Fe/H] £ — 1) show a linearly decreasing [O/Fe] with increasing metallicity 
(e.g. Edvardsson et al. 1993). The break near [Fe/H] ~ —1 was traditionally explained 
by the fact that the early heavy element production was done exclusively by SNe la, with 
the break occurring when the larger Fe production by SNe la kicks in (e.g. Matteucci & 
Greggio 1986). 

Recently, Kobayashi et al. (1998) performed chemical evolution calculations for both 
the double-degenerate scenario and for the single-degenerate scenario. For the latter they 
used two types of progenitor systems: one with a red giant companion and an orbital 
period of tens to hundreds of days, and the other with a near main sequence companion 
and a period of a few tenths of a day to a few days. 

They obtained for the double-degenerate scenario (for which they took a time delay 
of ~ 0.1-0.3 Gyr) a break at [Fe/H] ~ —2. For the single-degenerate scenario (with a 
delay caused by the main sequence lifetime of £ 1 Gyr; including metallicity effects), 
they obtained a break at [Fe/H] ~ —1. Kobayoshi et al. (1998) thus concluded that the 
double-degenerate scenario is inconsistent with Galactic chemical evolution results. 

Personally, I am not convinced by this apparent discrepancy, since Galactic chemical 
evolution calculations (and observations) are notoriously uncertain. In particular, the 
most recent Keck observations of oxygen in unevolved metal-poor stars show no break 
in the [O/Fe] vs. [Fe/H] relation. Rather, oxygen is enhanced relative to iron over three 
orders of magnitude in [Fe/H] in a robustly linear relation (Boesgaard et al. 1999). 
Consequently, apparent inconsistencies based on Galactic chemical evolution cannot be 
regarded at present as a fatal weakness of the double-degenerate scenario. 

5.1.2. SN la or accretion induced collapse? 

Potentially a more serious (and possibly even fatal) weakness of the double-degenerate 
scenario comes from the fact that some estimates and calculations indicate that the 
coalescence of two C-0 WDs may lead to an accretion-induced collapse rather than to 
a SN explosion (e.g. Mochkovitch & Livio 1990; Saio & Nomoto 1985, 1998; Kawai, Saio 
& Nomoto 1987; Timmes, Woosley k Taam 1994). 

The point is the following: once the lighter WD fills its Roche lobe, it is dissipated 
within a few orbital periods (Benz et al. 1990; Rasio & Shapiro 1995; Guerro 1994) 
and it forms a hot thick disk configuration around the more massive white dwarf. This 
disk is mainly rotationally supported and hence central carbon ignition does not take 
place immediately, but rather the subsequent evolution depends largely on the rate of 
angular momentum transport and removal, since they determine the accretion rate onto 
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the primary WD. As long as the accretion rate is higher than about M £ 2.7 x 10~ 6 M Q 
yr -1 , carbon is ignited off-center (at the core-disk boundary; this may happen during 
the merger itself; e.g. Segretain 1994). Under such conditions, the flame was found (in 
spherically symmetric calculations) to propagate all the way to the center within a few 
thousand years, thus burning the C-0 into an O-Ne-Mg mixture with no explosion (i.e. 
before carbon is centrally ignited; e.g. Saio & Nomoto 1998). Such configurations are 
expected to collapse (following electron captures on 24 Mg) to form neutron stars (Nomoto 
& Kondo 1991; Canal 1997). The main questions are then: 

(i) What accretion rates can be expected from the initial WD-thick disk configuration? 

(ii) May some aspects of the flame propagation be different given the fact that the 
real problem is three-dimensional while most of the existing calculations were performed 
using a spherically symmetric code? In particular, could the carbon burning be quenched 
before the transformation to O-Ne-Mg composition occurs? 

The answers to both of these questions involve uncertainties, however some possibilities 
are more likely than others. First, it appears very difficult to avoid high accretion rates. 
If the MHD turbulence that is expected to develop in accretion disks (e.g. Balbus & 
Hawley 1998) is operative, with a corresponding viscosity parameter of a ~ 0.01 (where 
the viscosity is given by v ~ ac s H, with H being a vertical scaleheight in the disk and 
c s the speed of sound; e.g. Balbus, Hawley & Stone 1996), then angular momentum can 
be removed in a matter of days! In such a case, even if the accretion rate is Eddington 
limited (at ~ 10~ 5 M /yr), off-center carbon ignition should still occur, with an eventual 
collapse rather than an explosion. Deviations from spherical symmetry can only hurt, 
since they may allow accretion to proceed at a super-Eddington rate. It is difficult to 
see why the dynamo-generated viscosity would be suppressed for the kind of shear and 
temperatures expected in the disk. 

Concerning the burning itself, recent attempts at multi-dimensional calculations of the 
flame propagation and a more detailed analysis of some of the processes involved (Garcia- 
Senz, Bravo & Serichol 1998; Bravo & Garcia-Senz 1999) indicate that if anything, 
accretion induced collapses are an even more likely outcome than previously thought. 
This is due to the effects of electron captures in Nuclear Statistical Equilibrium which 
tend to stabilize the thermonuclear flame, and to Coulomb corrections to the equation 
of state. The latter has the effect of reducing the flame velocities and the electronic and 
ionic pressures, all of which result in a reduction in the critical density which separates 
explosions from collapses. 

Finally, on the observational side there are also two points which argue at some level 
against WD mergers as SNe la progenitors. 

(i) Even if MHD viscosity could somehow be suppressed, and the disk surrounding the 
primary WD could cool down, so that angular momentum would be transported only via 
the viscosity of (partially) degenerate electrons, this would result in an accretion timescale 
of ~ 10 9 yrs (Mochkovitch & Livio 1990; Mochkovitch et al. 1997). The system prior to 
the explosion would have an absolute magnitude of My £10 (with much of the emission 
occurring in the UV). There is no evidence for the existence of some ~ 10 7 such objects 
in the Galaxy. 

(ii) The existence of planets around the pulsars PSR 1257+12 and PSR 1620-26 (Wol- 
szczan 1997; Backer 1993; Thorsett, Arzoumanian & Taylor 1993) could be taken to mean 
(this is a model dependent statement) that mergers tend to produce accretion induced 
collapses rather than SNe la. In one of the leading models for the formation of such 
planets (Podsiadlowski; Pringle & Rees 1991; Livio, Pringle & Saffer 1992), the planets 
form in the following sequence of events. The lighter WD is dissipated (upon Roche lobe 
overflow) to form a disk around the primary. As material from this disk is accreted, 
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matter at the outer edge of the disk has to absorb the angular momentum, thereby ex- 
panding the disk to a large radius. The planets form from this disk in the some way that 
they did in the solar system, while the central object collapses to form a neutron star. 

5.1.3. Overall assessment of the double-degenerate scenario 

It has now been observationally demonstrated that many double-degenerate systems 
exist. The general agreement between the distribution of the observed properties (e.g. 
orbital periods, masses) and those predicted by population synthesis calculations (Saffer, 
Livio & Yungclson 1998), suggests that the fact that no clear candidate (short period) 
system with a total mass exceeding the Chandrasekhar mass has been found yet, may 
merely reflect the insufficient size of the observational sample. Thus, there is very little 
doubt in my mind that statistics is not a serious problem. The most disturbing uncer- 
tainty is related to the outcome of the merger process itself. The discussion in §5.1.2 
suggests that collapse to a neutron star is more likely than a SN la (see also Mochovitch 
et al. 1997). 

5.2. The single-degenerate scenario 

The main strengths of the single degenerate scenario are: 

(1) A class of objects in which hydrogen is being transferred at such high rates that it 
burns steadily on the surface of the WD has been identified — the Supersoft X-ray Sources 
(Greiner, Hasinger & Kahabka 1991; van den Heuvel et al. 1992; Southwell et al. 1996; 
Kahabka and van den Heuvel 1997). If the accreted matter can indeed be retained, this 
provides a natural path to an increase in the WD mass towards the Chandrasekhar mass 
(e.g. Di Stefano & Rappaport 1994; Livio 1995, 1996a; Yungelson et al. 1996). 

(2) Other candidate progenitor systems are known to exist, like symbiotic systems 
(e.g. Munari & Renzini 1992; Kenyon et al. 1993; Hachisu, Kato & Nomoto 1999) and 
recurrent novae (Hachisu et al. 1999). 

(3) There have been claims that the single degenerate scenario fits better the results 
of Galactic chemical evolution (e.g. Kobayoski et al. 1998). However, as I have shown in 
§5.1.1, recent observations cast doubt on this assertion. Similarly, nucleosynthesis results 
show that in order to avoid unacceptably large ratios of 54 Cr/ 56 Fe and 50 Ti/ 56 Fe, the 
central density of the WD at the moment of thermonuclear runaway must be lower than 
~ 2 x 10 9 g cm~ 3 (Nomoto et al. 1997). Such low densities are realized for high accretion 
rates ( £ 1CP 7 M yr _1 ), which are typical for the Supersoft X-ray Sources. Nucle- 
osynthesis results suffer too, however, from considerable uncertainties (e.g. Nagataki, 
Hashimoto & Sato 1998). 

The main weaknesses of the single degenerate scenario are: 

(1) The upper limits on radio detection of hydrogen at 2 and 6 cm in SN 1986G, taken 
approximately one week before optical maximum (Eck et al. 1995), rule out a symbiotic 
system progenitor for this system with a wind mass loss rate of 1CT 7 < M w ~ 1CT 6 M Q 
yr _1 (Boffi & Branch 1995). This in itself is not fatal, since SN 1986G is somewhat 
peculiar (e.g. Branch and van den Bergh 1993), and the upper limit on the mass loss rate 
is at the high end of observed symbiotic winds. An even less stringent upper limit from 
x-ray and Ha observations exists for SN 1994D (Cumming et al. 1996). 

(2) There exists some uncertainty whether WDs can then reach the Chandrasekhar 
mass at all by the accretion of hydrogen (e.g. Cassisi, Iben & Tornambe 1998). Further- 
more, even if they can, the question of whether they can produce the required SNe la 
statistics is highly controversial (e.g. Yungelson et al. 1995, 1996; Yungelson & Livio 
1998; Hachisu, Kato & Nomoto 1999; Hachisu et al. 1999). 

I will now examine these weaknesses in some detail. 
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5.2.1. Observational detection of hydrogen 

Ultimately, the presence or total absence of hydrogen in SNe la will distinguish unam- 
biguously between single-degenerate and double-degenerate models. To date, hydrogen 
has not been convincingly detected in any SN la. It is interesting to note that narrow 
A6300, A6363 [01] lines were observed only in one SN la (SN 1937C; Minkowski 1939), 
but even in that case there was no hint of a narrow Ha line. Hachisu, Kato & Nomoto 
(1999) estimate in one of their models (which involves stripping of material from the red 
giant; see below) a density measure of M/vio ~ 10~ 8 M yr _1 (where v w is the wind ve- 
locity in units of 10 km s _1 ), while the most stringent radio upper limit existing currently 
(for SN 1986G) is M/v w ~ 10~ 7 M Q yr" 1 (Eck ct al. 1995; for SN 1994D Cumming 
et al. (1996) find from Ha an upper limit of M ~ 1.5 x 10~ 5 M Q yr _1 for a wind speed 
of 10 km s" 1 ; for SN 1992A Schlegel & Petre (1993) find from X-ray observations an 
upper limit of M/vio = (2 — 3) x 10~ 6 M Q yr -1 ). Thus, while it is impossible at present 
to rule out single-degenerate models on the basis of the apparent absence of hydrogen, 
the hope is that near future observations will be able to determine definitively whether 
this absence is real or if it merely represents the limitations of existing observations (an 
improvement by two orders of magnitude will give a definitive answer). 

5.2.2. Statistics 

Growing the WD to the Chandrasekhar mass is not easy. At accretion rates below 
~ 10~ 8 M /yr WDs undergo repeated nova outbursts (e.g. Prialnik & Kovctz 1995), 
in which the WDs lose more mass than they accrete between outbursts (e.g. Livio & 
Truran 1992). For accretion rates in the range 10 _8 -a few xlO~ 7 M Q /yr, while helium 
can accumulate, the WDs experience mass loss due to helium shell flashes and due to 
the common envelope phase which results from the engulfing of the secondary star in 
the expanding envelope (with mass loss occurring due to drag energy deposition). At 
accretion rates above a few xlO~ 7 M Q /yr, the WDs expand to red giant configurations 
and lose mass due to drag in the common envelope and due to winds (e.g. Cassisi et al. 
1998). The net result of this has been that population synthesis calculations which follow 
the evolution of all the binary systems in the Galaxy, tended until recently to conclude 
that single degenerate channels manage to bring WDs to the Chandrasekhar mass only 
at about 10% of the inferred SNe la frequency of 4 x 10~ 3 yr _1 (e.g. Yungclson ct al. 
1995, 1996; Yungelson & Livio 1998; Di Stefano et al. 1997; although see Li & van den 
Heuvel 1997). 

Very recently, a few serious attempts have been made to investigate whether the statis- 
tics could be improved by increasing the "phase space" for single degenerate scenarios, 
given the fact that population synthesis calculations involve many assumptions. These 
attempts resulted in three directions in which the phase space could be increased. 

(i) The accumulation efficiency of helium has been recalculated using OPAL opacities 
(Kato & Hachisu 1999). These authors concluded that helium can accumulate much more 
efficiently than found by Cassisi et al. (1998), mainly because the latter authors used 
relatively low WD masses (0.516 M Q and 0.8 M ) and old opacities in their calculations. 

(ii) Hachisu et al. (1999) claimed to have identified an evolutionary channel for single- 
degenerate systems previously overlooked in population synthesis calculations. In this 
channel, the C-0 WD is formed from a red giant with a helium core of 0.8-2.0 M Q 
(rather than from an asymptotic giant branch star with a C-0 core). The immediate 
progenitors in this case are expected to be either helium-rich Supersoft X-ray Sources or 
recurrent novae of the U Sco subclass. 
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(iii) It has been suggested that the inclusion of a few additional physical effects, can 
increase substantially the phase space of the symbiotic channel (Hashisu, Kato & Nomoto 
1999). These new effects included: 

(1) The WD loses much of the transferred mass in a massive wind. This has the 
effect that the mass transfer process is stabilized for a wider range of mass ratios, up to 
Qmax = m 2 /mi =1.15 instead of g max = 0.79 without the massive wind. 

(2) It has been suggested that the wind from the WD strips the outer layers of the 
red giant at a high rate. This increases the allowed mass ratios (for stability) even above 
1.15, essentially indefinitely. 

(3) It has been suggested that at large separations (up to ~ 30, 000 R Q ) the wind from 
the red giant acts like a common envelope to reduce the separation, thus allowing much 
wider initial separations to result in interaction. 

There are many uncertainties associated with all of these attempts to increase the phase 
space. For example, the efficiency of mass stripping from the giant by the wind from 
the WD may be much smaller than assumed by Hachisu et al. (1999), for the following 
reasons. At high accretion rates, much of the mass loss from the WD may be in the form 
of an outflow or a collimated jet, perpendicular to the accretion disk rather than in the 
direction of the giant. Evidence that this is the case is provided by the jet satellite lines 
to He II 4686, H/3 and Ha observed in the Supersoft X-ray Source RX J0513.9-6951 
(Southwell et al. 1996). These jet lines are very similar to those seen in the prototypical 
jet source SS 433 (e.g. Vermeulen et al. 1992). Furthermore, even if some of the WD 
wind hits the surface of the giant, it is not clear how efficient it would be in stripping 
mass, since the rate of energy deposition per unit area by the wind is smaller by two 
orders of magnitude that the giant's own intrinsic flux. 

Similarly, the efficiency of helium accumulation is still highly uncertain, as the dif- 
ferences between the results of Kato & Hachisu (1999) and Cassisi et al. (1998) have 
shown. 

Finally, all the new suggestions for the increase in phase space rely very heavily on 
the results of the wind solutions of Kato (1990; 1991), which involve a treatment of the 
radiation not nearly as sophisticated as that of more state of the art radiative transfer 
codes (e.g. Hauschildt et al. 1995, 1996). 

5.2.3. Overall assessment of the single-degenerate scenario 

The above discussion suggests that probably not all the scenarios for increasing the 
"phase space" of the single-degenerate channels work (if they did we might have had the 
opposite problem of too high a frequency of SNc la!). However, these attempts serve 
to demonstrate that the input physics to population synthesis codes still involves many 
uncertainties. My feeling is therefore that given the many potential channels leading to 
SNe la, statistics should not be regarded as a serious problem. 

Single-degenerate scenarios therefore appear quite promising, since unlike the situation 
a decade ago, a class of objects in which the WDs accrete hydrogen steadily (the Supersoft 
X-Ray Sources) has actually been identified. The main problem with single-degenerate 
scenarios remains the non-dctcction of hydrogen so far. While a difficult observational 
problem (see §6), the establishment of the presence or absence of hydrogen in SNe la 
should become a first priority for SNe observers. 

5.3. What if....? 

Given the fact that there are still uncertainties involved in identifying the SNe la progeni- 
tors, and that WD mergers and some form of off-center helium ignitions almost certainly 
occur, it is instructive to pose a few "what if" questions. For example: What if WD 
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mergers with a total mass exceeding Chandrasekhar do not produce SNe la, what do they 
produce then? The answer in this case will have to be that they almost certainly pro- 
duce either neutron stars via accretion induced collapses, or single WDs, if the merger is 
accompanied by extensive mass loss from the system. 

What if off-center helium ignitions do not produce SNe la? In this case, if an explosive 
event indeed ensues, a population of "super novae" (with ^0.15 M of 56 Ni and He) is 
yet to be detected (maybe SN 1885A in M31 was such an event?). What if off-center 
helium ignitions do produce SNe la? What comes out of the systems with M WD £ 1 M Q , 
which should be even brighter? It is difficult to believe that the latter are represented by 
the very few bright objects like SN 1991T. Thus, we see that off-center helium ignitions 
seem to present an observational problem both if they do and if they do not produce 
SNe la. To me this suggests that the physics of these events is not well understood (for 
example, maybe off-center helium ignition fails to ignite the C-0 core after all). 

6. How can we hope to identify the progenitors? 

There are several ways in which observations of nearby supernovae could solve the 
mystery of SNe la progenitors: 

(1) A combination of early high resolution optical spectroscopy, x-ray observations and 
radio observations can both provide limits on M/v from the progenitors and potentially 
detect the presence of circumstellar hydrogen (if it exists). 

For example, narrow HI in emission or absorption could be detected either very early, 
or shortly after the ejecta become optically thin (~ 100 days). The latter is true because 
the SN ejecta probably engulfs the companion at early times (e.g. Chugai 1986; Livne, 
Tuchman & Wheeler 1992). The interaction of the ejecta with the circumstellar medium 
can be observed either in the radio (e.g. Boffi & Branch 1995) or in x-rays (e.g. Schlegel 
1995). The collision of the ejecta (with circumstellar matter) can also set up a forward 
and a reverse shock (e.g. Chevalier 1984; Fransson, Lundqvist & Chevalier 1996), and 
radiation from the latter can ionize the wind and produce Ha emission (e.g. Cumming 
et al. 1996). 

(2) Early observations of the gamma-ray light curve (or gamma-ray line profiles) could 
distinguish between carbon ignitors and sub-Chandrasekhar helium ignitor models (see 
§4.2.2) since the latter can be expected to result in a quicker rise of the gamma-ray 
light curve due to the presence of 56 Ni in the outer layers (and different gamma-ray line 
profiles; because of the high velocity 56 Ni). 

Observations of very distant supernovae (at z ~ 3-4) with the Next Generation Space 
Telescope (NGST) can also help (e.g. Yungelson & Livio 1999). For example, the pro- 
genitors can be identified from the observed frequency of SNe la as a function of redshift 
(e.g. Yungelson & Livio 1998, 1999; Ruiz-Lapuente & Canal 1998; Madau, Delia Valle 
& Panagia 1998; and see §8), since different progenitor models produce different redshift 
distributions. Personally, I think that it would be absolutely pathetic to have to resort 
to this possibility. Rather, one would like to identify the progenitors independently, and 
to use the observations of supernovae at high z to constrain models of cosmic evolution 
of rates, luminosity, and input into galaxies. 

7. Could we be fooled? 

One of the key questions that result from the fact we do not know with certainty 
which systems are the progenitors of SNe la is clearly: is it possible that SNe la at higher 
redshifts are systematically dimmer than their low-redshift counterparts? In this respect 
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it is important to remember that a systematic decrease in the brightness by ~ 0.25 
magnitudes is sufficient to explain away the need for a cosmological constant. In a recent 
work, Yungelson & Livio (1999) calculated the expected ratio of the rate of SNe la to 
SNe from massive stars (Types II, la, Ic) as a function of redshift for several progenitor 
models. They showed that if different progenitor systems can contribute to the total 
SNe la rate (e.g. double-degenerates at the Chandresekhar mass and single-degenerates 
with subgiant donors at sub-Chandrasekhar masses), then it is possible in principle that 
a different progenitor class will start to dominate at z ~ 1. However, such a transition 
is highly unlikely, because: (i) in some models the transition has the opposite effect to 
the observed one (e.g. double degenerates which may be expected to be brighter than 
sub-Chandrasekhar ELDs dominate at the higher redshifts). (ii) If the contribution from 
physically different channels was indeed significant, one would have expected to observe 
this division more clearly also in the local sample, which is not the case (~ 90% of SNe la 
are "normals"). Consequently, I do not believe that the observed universal acceleration 
is an artifact of the observed sample being dominated by different progenitor classes. 



8. Tentative conclusions and observational tests 

On the basis of the analysis and discussion in the present work, the following tentative 
conclusions can be drawn: 

(1) SNe la are almost certainly thermonuclear disruptions of mass accreting C-0 white 
dwarfs . 

(2) It is very likely that the explosion occurs at the Chandrasekhar mass, as carbon 
is ignited at the WD center. Off-center ignition of helium at sub-Chandrasekhar masses 
may still be responsible for a subset of the SNe la which are subluminous, but this is less 
clear. 

(3) The immediate progenitor systems are still not known with certainty. From the 
discussion in §5 (see in particular §5.1.3 and 5.2.3) however, I conclude that presently 
single degenerate scenarios look more promising, with hydrogen or helium rich material 
being transferred from a subgiant or giant companion (systems like Supersoft X-Ray 
Sources and Symbiotics). 

(4) Definitive answers concerning the nature of the progenitors can be obtained from 
observations taken as early as possible in: x-rays, radio, and high resolution optical 
spectroscopy. The establishment of the presence or absence of hydrogen in SNe la should 
be regarded as an extremely high priority goal for supernovae observers. If hydrogen will 
not be detected at interesting limits (corresponding to M /v w <~ 10~ 8 M yr _1 ), this 
will point clearly towards the double-degenerate scenario. 

(5) Observations of SNe la at high redshifts can help to test particular ingredients of 
the models which are directly related to the nature of the progenitors. For example, most 
of the models aiming at improving the statistics of the single-degenerate scenarios rely 
on a strong wind from the accreting WD. These models thus predict an "inhibition" of 
SNe la in low-metallicity environments, and in particular a significant decrease in the rate 
of SNe la at z <~ 1-2 (Kobayashi et al. 1998). At present, the detection of a very likely 
SN la at redshift z = 1.32 (SN 1997ff) in the Hubble Deep Field appears inconsistent 
with this prediction (Gilliland, Nugent & Philips 1999), but more observations will be 
required to give a more definitive answer. 

This research has been supported in part by NASA Grant NAG5-6857. 
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